Electronic transport properties through thiophenes on switchable domains 
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The electronic transport of electrons and holes through stacks of Q:,ciJ-dicyano-/3,/3'-dibutyl- 
quaterthiophene (DCNDBQT) as part of a novel organic ferroic field-effect transistor (OFFET) 
is investigated. The novel application of a ferroelectric instead of a dielectric substrate provides 
the possibility to switch bit-wise the ferroelectric domains and to employ the polarization of these 
domains as a gate field in an organic semiconductor. A device containing very thin DCNDBQT 
films of around 20 nm thickness is intended to be suitable for logical as well as optical applications. 
We investigate the device properties with the help of a phenomenological model called multilayer or- 
ganic light-emitting diodes (MOLED), which was extended to transverse fields. The results showed, 
that space charge and image charge effects play a crucial role in these organic devices. 

PACS numbers: 72.80.Le, 77.84.-s, 81.07.Pr, 87.15.hj, 



I. INTRODUCTION 

The investigation of organic semiconductors, which 
consist of special organic molecules or polymers, has 
become more and more attractive in recent year^^Il 
From the experimental point of view, the main fo- 
cus was directed to the fabrication and characteriza- 
tion of suitable organic devices such as organic light- 
emitting diodes (OL EDs) and organic field-effect tran- 
sistors (OFETs)EEsl. Theoretical effort was linked to the 
description of the physical background, which includes 
the investigation of the underlying transport mecha- 
nisms. A mobile charge carrier can be transported 
ballistically along a polymer or it jumps between suit- 
ably stac ked molecule s, which is well-known as hopping 
transporii ^ l ^ l ^ l ^ l ^^ l -'^^ l ^^ l In addition, the characterization 
of the orga nic-ele ctrode interface is essential for potential 
application^i^IIl] 

We present theoretical investigations of a novel or- 
ganic device, which can be formally denoted as an organic 
ferroic field effect transistor (OFFET). Contrary to the 
dielectric substrates generally employed in conventional 
OFET setups (see upper panel of figure [l]), OFFET de- 
vices are fabricated with ferroelectric materials such as 
BaTiOs or Pb(Zr,Ti)03. Below a material-dependent 
critical temperature T^, the ferroelectric material spon- 
taneously polarizes and forms ferroelectric domains with 
different directions of polarization. In the OFFET this 
induced polarization field serves as gate field. Due to the 
different operation principle, this can be seen as an alter- 
native to the well-established high perfor mance ferroelec- 
tric field-effect transistor (FcFET) devicePEsE^EHUnilol^ 

The ideal ferroelectric substrate would form domains 
with two preferable directions of polarizibility. Nowa- 
days the orientation of these ferroelectric domains can 
be resolved for instance with the help of near-field 
microscopj^. The ideal ferroelectric substrate would 



form domains with two preferable directions of polariz- 
ability, in-plane and orthonormal. An in-plane polar- 
ization is expected to have only minor infiuence on the 
organic medium, whereas an polarization field Ftr in or- 
thonormal direction is supposed to have maximum infiu- 
ence. Consequently, the organic material must be po- 
larizable in this orthonormal direction. Therefore, it is 
essential for potential applications to understand the in- 
fluence of the ferroelectric domains on the organic mate- 
rial. Thus, we will study the polarizability of the organic 
medium. 




FIG. 1: [left] Enlarged slice of the thiophene film represent- 
ing its morphologj^^. [right] General setup of an OFET device 
(upper panel) and the modified version (lower panel), where 
the dielectric substrate is replaced by a ferroelectric material 
such as PZT. The ferroelectric domain induces a polarization 
field which influences the organic material. Grey and black 
arrows at the bottom represent different directions of polar- 
ization. 

The employed organic material will be a,a;-dicyano- 
/3,/?'-dibutylquaterthiophene (DCNDBQT). This aro- 
matic molecule, belonging to the group of the oligoth- 
iophenes, can form well-ordered stacks with a good tt- 
orbital overlap. This high degree of ordering can favor 
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high carrier mobihties and additionally decrease disorder 
effects, which are typically observed in disordered organic 
medicPl. The DCNDBQT films show thermal stability 
after preparation with an energy gap Eg in the range of 
semiconducting materials {Eg < 3eV)^. 

Furthermore, we will present the investigation of the 
electronic transport in such an OFFET device. For 
this purpose, the device is simulated with the help of a 
phenomenological numerical model called multilayer or- 
ganic light-emitting diode (MOLED )P2l2ll and its exten- 
sion to transverse fields. While transport properties can 
be found elsewher^^i], we will now focus on electric prop- 
erties influenced by the ferroelectric gate. In addition, 
the simulation-relevant parameters are determined either 
by experimental, analytical or computational methods. 

Another crucial point regards the carrier transport in 
organic systems, which can be very complex due to dy- 
namic polarization effects, either in the organic medium 
itself and additionally in the dielectric substrat^^. Due 
to the low carrier densities in our system, these intrinsic 
polarization effects will play only a minor role. Addi- 
tional experiments are necessary to clearify this assump- 
tion. 

Generally, the magnitude of the gate voltage is con- 
stant and depends on the employed ferroelectric material. 
Thus, the gate voltage in our simulations is a simula- 
tion parameter, which has later to be identified with the 
intrinsic, material-dependent surface polarization. The 
material-intrinsic electric field can be roughly estimated 



by employing Ft, 
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under the assumption of a ho- 



mogeneously charged ferroelectric surface. Here, a repre- 
sents the surface charge, which lies typically in the range 
of 30^ for materials like Pb(Ti,Zr)03 and LiNbOs, 
while er — 2.8 denotes the relative dielectric constant 
of the molecular crystal of thiophen^^H The estimated 
polarization field amounts to Ftr ~ directly at the 
interface between substrate and organic. Furthermore, 
a decrease of this electric field is expected further inside 
the organic medium due to screening effects. This might 
limit the organic film thickness, particularly in logical de- 
vices, where switching between two states is important. 

Ferroganic devices have a great potential for optical 
as well as logical applications. While optical ferroganic 
devices are interesting due to the possibility of band-gap 
engineering (see chapter 2 for detailed explanation), OF- 
FET structures take benefit of the domain stability in a 
working device. This advantage is given by the material- 
specific gate field which persists even after a power out- 
age. 

Section |TT] presents the influence of the polarization 
field on the molecule. Section |III| gives some informa- 
tion about the alorithm employed in the MOLED model 
and the determination of the simulation parameters. Sec- 
tions IIVI and |V] show the results of the simulation and a 
conclusion, respectively. 



II. MOLECULE POLARIZATION 

The modeling of the OFFET requires a detailed under- 
standing of how the molecule DCNDBQT will be influ- 
enced by the application of a transverse electric field Ftr, 
which is generated by the ferroelectric material. There- 
fore, the field-induced variation of the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) is of particular interest be- 
cause these orbitals correspond to the valence and con- 
duction band of inorganic semiconductors, respectively. 





FIG. 2: DFTB-calculated HOMO deformation of DCNDBQT 
by the influence of a transverse electric field in positive x di- 
rection. The different colors denotes the different signs (+/-) 
of the electron wave function, (a): DCNDBQT with coordi- 
nate system, (b)-(d): Increase of the field strength from zero 
field, over Fi = 2 ^ up to Fa = 6 
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FIG. 3: Variation of the HOMO/LUMO energy levels of 
quaterthiophene under the application of a transverse field 
Ftr- The band gap decreases linearly (region I) with in- 
creasing strength of the electric field. The non-linear regime 
(region II) is reached at fields higher than 6 when the 
HOMO-LUMO compositions change. 

For this purpose, the influence of the transverse 
electric field on the electronic properties of oligothio- 
phenes has been studied using the self-consistent charge 
densit y-functi onal based tight binding method SCC- 
DFT^2Z128i2l!. Figure H] illustrates, that the 7r-electron 
system is most strongly deformed along the molecule with 
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increasing i^trll^;, hence polarizing the molecule towards 
one short edge. 

Figure |3] shows the resulting variation of the HOMO 
and LUMO energy levels of DCNDBQT under an applied 
transverse electric field Ft^-Ha;. This field Ftr induces a 
compensating polarization in the molecule, which corre- 
sponds to a mixing of HOMO and LUMO levels. This 
leads to a linear decrease of the undistorted HOMO- 
LUMO gap, in detail from 1.8 eV at zero field to 0.6 
eV at a high external field of 6 — (region I). Conse- 
quently, this molecule is suitable for devices which take 
advantage of a band-gap engineering. 

Figure |3] also shows a transition to non- linear behav- 
ior at fields higher than 6 ^ (region H). For our sim- 
ulations, we will stay in the linear regime by employing 
fields not higher than Ffr = 4 ^ . 

III. NUMERICAL DETAILS 

We employ the phenomenological model MOLED, in- 
troduced by Houili et al.^'', to investigate the charge car- 
rier profile in our OFFET device. A main advantage of 
this model is its fiexibility due to the full parametriza- 
tion of the device. These parameters comprise, for in- 
stance, the electrode work functions, molecular energy 
levels, temperature and applied electric fields. 

In this section, we will present first the general con- 
cepts and boundary conditions of MOLED and second, 
the determination of the most important simulation pa- 
rameters. 



A. General concepts and boundary conditions of 
MOLED 

MOLED is based on an effective 1-dimensional model, 
which can be applied, if the system fulfills three key as- 
sumptions, (i) The electrodes are planar and parallel to 
each other, (ii) The organic channel material between 
the electrodes is homogeneously distributed and parallel 
to the electrodes, (iii) The distance between the elec- 
trodes should be small compared to the other two spatial 
dimensions. 

Experimentally, the first point can be achieved by 
photolithographic fabrication of the electrodes. Fur- 
thermore, the AFM image in ref. El depicts the 
high degree of order of the DCNDBQT molecules 
in the film, which is promoted by a template layer 
consisting of 5-cyano-2-(butyl(4-phosphonic acid))-3- 
butylthiophene (CNBTPA). Therefore the second as- 
sumption is also fulfilled. The authors of ref. ^5 deter- 
mined the molecule-molecule distance to 3.5 A, confirmed 
by a self-consistent charge density- functional based tight- 
binding calculation of the molecule-molecule equilibrium 
distanc^^. The third assumption of MOLED is not 
fulfilled because the DCNDBQT film is only a couple 
of monolayers thick. However, the AFM data depicts 



an array of well-ordered chains of 7r-stacked DCNDBQT 
molecules, hence they indicate that the transport occurs 
through parallel quasi-lD channels with little cross-talk. 

Although this 1-dimensional approach is very straight- 
forward it can lead to results different from the simulation 
of real 3-dimensional systems. MOLED assumes a homo- 
geneous charge distribution within the molecule. Strictly 
speaking, this is not the case, especially close to the elec- 
trodes. Moreover, charge carriers may also be more lo- 
calized at a certain region on the molecule due to pola- 
ronic effects, consequently leading to inhomogeneities of 
the charge distribution. Thus, future investigations are 
necessary to clarify this situation in our system. 

For the modeling of an organic semiconductor, a de- 
tailed knowledge of its intrinsic conductance is required. 
In MOLED, three key mechanisms are included, which 
capture the physics of an organic semiconductor: First, 
the charge carrier tunneling from the electrodes into 
the channel, second the charge carrier transport via a 
phonon-assisted hopping process and third, charge car- 
rier recombination and emission of photons. 
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FIG. 4: Sketch of MOLED setup. Molecules are situated 
along the simulation axis between the cathode and anode. 
Electrons are injected at the cathode and holes at the anode. 
The charge carriers move along the molecule chain via a ther- 
mally assisted hopping process. The HOMO and LUMO of 
the molecule correspond to the conduction and valence bands, 
respectively. Recombination occurs, when an electron and 
hole are within a characteristic distance. 

Figure |4] gives an overview of the simulation setup and 
depicts these three mechanisms. The MOLED simulation 
axis is ahgned perpendicularly to the electrodes and the 
hopping sites, which correspond to single molecules, are 
discrete in space. Thus, the molecules can be seen as flat- 
shaped objects on a 1-dimensional chain. Furthermore, 
MOLED is a two- niveau energy level approach, consider- 
ing the HOMO and LUMO energy levels of the molecules 
for hole and electron transport, respectively. 

The three main mechanisms of our simulated organic 
semiconductor are based on these molecular energy lev- 
els, which can be modified, for instance by Coulomb in- 
teractions, an external electric field, image forces and 
space-charge effects. 

Consequently, the electronic states of an organic 
molecule can be characterized by: 
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Parameter 



value 



Em = Eo,,n + Eif{x„i) + EcsjXm) + ED{Xm\ (1) 
reference variation 

with m ~ 1, ...,N . Equation[l]detcrmincs the value of the 
molecular energy level of a molecule at a given position 
Xm- The first term £^o,m is the reference term, which cor- 
responds to the bare energy levels of the organic molecule 
without any external influence factors. The HOMO and 
the LUMO energies are used as reference values for holes 
p and for electrons e, with site-specific concentrations Pm 
and I respectively. These bare energy levels can be de- 
termined by experimental methodd^ or with the help of 
DFT calculationsP. The remaining terms of eq. [l] ac- 
count for the infiuence of the local environment of the 
molecule and are denoted as the variation term. In de- 
tail, the first term Eip{xm) takes into account the image 
force potential. The next term, Ecsi^m) represents the 
Coulomb shift of the energy levels caused by the inter- 
action of the homogeneously charged sheets with each 
other. Er){xm,) is a correction term, which models an in- 
homogeneous charge distribution close to the electrodes 
instead of the general homogeneous distribution inside 
the device. For more details, the reader is referred to ref. 
[251 and to the references therein. 



"Distance between molecules 

"Distance molecule <-> electrode 

"HOMO 

"LUMO 

"^AE (HOMO) 

=AE (LUMO) 

"Work function of electrodes (Au) 

"DCNDBQT height 

Discretization length r 

Simulation temperature 

"Relative permittivity of thiophene 

"Hole mobility fig 

''"'^Tunneling parameter 7^'^ ■ <?2'fl 

Tunneling prefactor A" 

Tunneling prefactor A'' 

Attenuation factor a" 

Attenuation factor 

Mobility law 

''Poole-Frenkel factor _fo 

Convergence parameter 

Voltage interval [Vmin;Vmax;Vsteps] 



3.5 A 

3.0 A 
-5.58 eV 
-3.33 eV 

0.096 

GV/m 

0.113 

GV/m 

5.1 eV 
15 A 

8.74 A 
300 K 
2.8 

1.6 • 10~5 ^ 

Vs 

0.73 electrons 
3.3 A 
2 A 
0.7 
0.7 

Poole-Frenkel 
103416 f 
l.OE-12 

[-3;-12;50] 



TABLE I: Overview of the employed MOLED parameters. 
The parameters are determined with experimental a, analyti- 
cal b or computational c methods. The remaining parameters 
are taken from the reference modePl. 



B. Determination of simulation- relevant 
parameters 

A crucial part of the model are the system-relevant 
OFFET parameters, for which realistic values have to be 
determined. Table |l] gives an overview of the employed 
parameter set. 

The parameters of category a are taken from the ex- 
perimental paper of ref. [5] Hopping is assumed between 
neighboring sites in reasonable agreement with the tight- 
binding approach of Morawetz et al.Hl. The underlying 
transport process follows the Poole-Frenkel mobility law 
with a calculated Fq as described in ref. |23] and [SO] The 
zero-field mobilities were taken from ref [51 The remain- 
ing section will present a detailed discussion about the 
determination of the tunneling parameter. It is calcu- 
lated with the help of the density of states (DOS) gL,R 
of the right and left electrode and the tunneling factor 

The DOS is difficult to access experimentally. A suit- 
able access to this quantity is provided in the framework 
of DFT, here in the plane- wave pseudopotential approach 
of the program package abinitISi]. With such calcula- 
tions it is possible to determine the DOS of the metal 
surface and thus estimate the contribution of mobile elec- 
trons to the transport from the electrode to the molecules 
at the interface. The upper panel of figure [5] depicts the 
shape of the calculated DOS per gold atom of the unre- 
constructed Au(lll) surface (dotted curve) , which is the 
preferred surface in nanoscale structures^Sl. 



The most important contribution arises from the elec- 
trons which have an energy close to the Fermi level Ep- 
To include temperature effects one can convolute the cal- 
culated DOS per atom D{E) with the Fermi-Dirac func- 
tion at the expected device temperature of 300 K. The 
resulting curve DrediEjis also shown as solid line in the 
upper panel of figure Kn The value g^\, which repre- 
sents the fraction of electrons provided by each surface 
gold atom, can now be obtained by integrating the func- 
tion Dred{E) over the region modified by the Fermi-Dirac 
function, consequently 

/•OO 

91%^ / Dred{E)dE. (2) 

To estimate the whole electron contribution to the tun- 
neling, one has to determine the number of those atoms, 
which are accessible for the transport over the electrode- 
molecule surface. This can be achieved in a first order ap- 
proximation by a simple geometrical estimate. The lower 
panel of figure [5] depicts the geometric situation at the 
interface. If one considers both, the lattice parameter of 
gold at the surface and the area covered by the molecule, 
one can determine the number of gold atoms Nau in the 
vicinity of the molecule, which leads to Nau = 11. Thus, 
each gold surface atom, which is situated directly below 
the TT-orbital system of DCNDBQT, can contribute to 
the tunneling process. MOLED simulations with differ- 
ent values of Nau show, however, that small variations in 
the total number of participating gold atoms have only 
a minor influence on the injection term, which indicates. 
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FIG. 5: [top] Density of states (DOS) per surface atom of an 
unreconstructed Au(lll) surface. Energies are shifted such 
that zero energy corresponds to the Fermi level. The dot- 
ted curve represents the calculated DOS. The solid curve is 
a convolution of the calculated DOS employing the Fermi- 
Dirac function at a temperature of 300K. This accounts for 
thermally activated electrons, [bottom] Geometrical repre- 
sentation of the electrode - molecule interface and participat- 
ing surface gold atoms to the tunneling process (highlighted 
grey). 



that the model of the electrode-molecule interface is in- 
sensitive against small geometric variations. The total 
value of gL,R can be determined by 



(3) 



The second parameter jl,r is also difficult to determine 
experimentally, because the employed molecule DCND- 
BQT has a negligible tunneling resistancd^. Instead, a 
value can be deduced from an analytical solution, where 
the Schottky barrier is calculated for a mean source-drain 
field of F = 0.65^. 

Here, the electrode is considered as a confining poten- 
tial well, which restricts the motion of the mobile elec- 
trons. During the injection process, the electron has to 
overcome a characteristic potential V to tunnel directly 
from the metal to the molecule. In this approach, the 
transmission coefficient T, which denotes the probability 



of the electron transmission, can be calculated analyti- 
cally as 



T = 



. 2 

sinh (k) 



AE{V - E) 



(4) 



where E is the energy of the electron wave, V the poten- 
tial barrier between electrode and molecule, and 



byf2m{V ~ E) 



(5) 



m is the mass of the electron and b the tunneling distance. 
From the assumption that k, fulfills k ^ 1, it follows for 
the transmission coefficient T 



WE{V -E) 

y2 



exp(— 2k). 



(6) 



The prefactor of the exponential term can be identified 
with the simulation parameter "fL,R (see for instance eq. 
5 of ref[23|, consequently 



1L,R 



IQE{V -E) 



(7) 



The remaining task is to determine the potential barrier 
V , which is a function of the undistorted barrier <i>B of 
the molecule-electrode interface, the image force (second 
term) and the field contribution (third term in eq |8| . 
Consequently, 



V= fi^B- 



qFx)dx 



(8) 



had to be integrated over the electrode - molecule dis- 
tance of 6 = 3 A. This electrode - molecule distance is 
in agreement with literature data^. In the physical pic- 
ture, the electron energy E for tunneling through the 
barrier can be approximated by the Fermi energy Ep of 
the metal electrode. 



IV. RESULTS 

In this section we present the results on the OFFET 
performance for a device with an electrode-electrode dis- 
tance of 1000 A at different gate voltages of Vg = {0; 1; 2} 
volts. These voltage values represent an effective polar- 
ization field of Ftr — {0;0.5; 1}^^^ in our ID thiophene 
stack in accordance to the experimental electric field^. 

The device distance was chosen to model a realistic 
extend of the ferroelectric domain of the gat^^. Thus, 
the switching of the ferroelectric substrate will polarize 
the whole area between the leads. We will investigate 
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FIG. 6: HOMO and LUMO energy levels for different gate 
voltages Vg in absence (upper panel) and presence (lower 
panel) of an applied source-drain voltage of Vsd ~ 10 volts. 
The Fermi level of the gold electrode is also shown. 



quantitatively the influence of the polarization on the 
charge carrier density at constant source-drain field. 

Figure [6] depicts the HOMO and LUMO energy levels 
without (upper panel) and with (lower panel) an applied 
source-drain voltage of Vsd ~ 10 volts. 

Although the energy levels are localized molecular 
states and do not form a band, a continous representa- 
tion was chosen, because the device is much larger than 
the intermolecular spacing. Both panels of figure [6] shows 
also the position of the Fermi level of the gold electrodes. 
Obviously, the injection barrier Eb for holes, which en- 
ters from the left side into the device, is much lower in 
comparison to the electron injection barrier on the cath- 
ode side. We estimated Eb for holes as E^ « 0.3 eV and 
for electrons as i?^ « 1.8 eV, respectively. Consequently, 
the holes can be considered as the majority charge carrier 
class, in accordance with experiment^. 

An applied gate voltage Vg will decrease the HOMO- 
LUMO energy gap, as indicated in the upper panel of 
figure [6] We expect, that there might occur an increase 
of the hole carrier density inside the device due to the 



lowering of the energy barrier for injection. Indeed, a 
closer look at the first molecules near the anode confirm 
an increase of the hole carrier density with increasing 
gate voltage Vg, as depicted in the upper right panel of 
figure [7| 
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FIG. 7: Hole density (majority charge carrier class) depending 
on the device position. [Top] A magnification of the anode- 
organic interface. [Bottom] The hole density of the whole 
device. 

However, the lower panel of figure [7] shows a decrease 
of the hole density further inside the device for a gate 
voltage of = 2 volts. This can be explained by the im- 
age force contribution, which arises at the metal-organic 
interface. The accumulated holes at the interface intro- 
duce image charges within the electrode, which act as an 
effective Coulomb blockade. Consequently, a counteract- 
ing field against the source-drain bias is created, which is 
shown in figure [8] This effect plays a minor role for gate 
voltages below or equal to = 1 volt. 

Additionally, this effect becomes pronounced due to 
the low hole carrier mobility of fj,o ^ 10~^^, by which 
the charge carriers are not transported away sufficiently 
fast from the injecting electrode. 

Figure [Ojdepicts the HOMO energy levels of the first 26 
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FIG. 8: [left] Local electric field in dependence of the applied 
gate voltage. 
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FIG. 9: [right] Variation of the Schottky barrier due to the 
change of the image force contribution. 



V. SUMMARY &: OUTLOOK 

The electronic transport of electrons and holes through 
stacks of a,a;-dicyano-/3,/3'-dibutyl- quaterthiophene as 
part of a novel organic ferroic field-effect transistor was 
presented. 

We investigated the influence of the ferroelectric po- 
larization of the substrate on the electronic transport 
through the organic layer with the help of the phe- 
nonicnological model MOLED, which treats the hopping 
transport of electrons and holes in such films by a mas- 
ter equation approach. To accomplish that, an extension 
of MOLED to include the influence of transverse electric 
fields was essential. The effect of the transverse field can 
be easily modeled by a linear decrease of the HOMO- 
LUMO gap for moderate field strengths, which we have 
calculated using the SCC-DFTB method. 

A subtask of this study was the determination of the 
simulation-relevant parameters by experimental, analyt- 
ical and computational methods. We presented detailed 
approaches for the determination of two essential simu- 
lation parameters, the density of states at the gold elec- 
trode and the tunnehng prefactor for injection into the 
organic medium. Both parameters are hardly accessible 
by experimental methods for this special kind of device. 

After determining the MOLED parameters, we quali- 
tatively studied the performance of a fully polarized OF- 
FET device. We investigated the device behavior for 
different gate voltages and showed the capability of the 
model for organic device analysis. 

Further investigations will focus on the comparison 
with prototypical devices and a detailed understanding 
of the underlying physics. 



molecules near the anode-organic interface. Obviously, 
the change in the image force contribution influences the 
shape and position of the Schottky barrier. This Schot- 
tky barrier is slightly elevated, which increases the ef- 
fective energy barrier for charge carrier injection, conse- 
quently leading to a decrease of the hole carrier density 
further inside the device. 
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